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1  Introduction 


Thermographic  mapping  of  microwave  fields  is  a  method  of  imaging 
the  electromagnetic  field  in  a  region  of  space  by  viewing  the  heating  it 
produces  when  interacting  with  a  material.  Many  people  are  familiar  with 
the  effects  of  microwaves  on  materials  through  the  wide-spread  use  of 
microwave  ovens.  As  the  microwave  field  interacts  with  a  lossy  material, 
some  of  the  electromagnetic  energy  may  be  converted  into  heat.  This  heat 
can  be  viewed  in  several  ways  [1],  the  most  useful  technique  developed  to 
date  involves  an  infrared  camera. 

This  method  of  measuring  microwave  fields  has  two  major  advantages 
over  conventional  methods  involving  metallic  probes.  First,  the  field-of- 
view  of  the  camera  is  usually  digitized  into  at  least  128  by  128  independent 
measurements  of  the  microwave  field,  yielding  more  information  in  a 
shorter  time.  Second,  the  material  which  heats  in  the  microwave  field  is  less 
perturbing  to  that  field  than  the  introduction  of  a  metallic  probe,  yielding 
more  accurate  measurements.  The  primary  limitation  up  to  this  time  has 
been  an  inability  to  associate  a  microwave  field  strength  with  an  observed 
temperature. 

The  objective  of  this  work  is  to  develop  appropriate  theoretical  and 
experimental  practices  to  allow  for  the  design  and  manufacture  of  an 
electromagnetic  absorbing  screen  which  will  allow  an  infrared  image  to  be 
converted  into  a  microwave  field  strength. 

An  exhaustive  development  of  the  pertinent  electromagnetic  theory  is 
given  first.  This  is  then  applied  to  the  problem  of  optimizing  the  dielectric 
constant,  conductivity  and  thickness  of  a  planar  absorber  in  order  to 
maximize  the  power  absorption  under  plane  wave  illumination.  A 
particular  variety  of  absorber,  a  carbon  or  ferrite  loaded  epoxy,  is  then 
manufactured.  A  method  is  developed  to  measure  the  complex  permittivity 
and  permeability  of  the  constructed  absorber  at  microwave  frequencies.  The 
end  result  of  this  work  is  an  absorber  with  known  electrical  characteristics. 

This  work  addresses  the  electrical  characteristics  of  an  absorbing 
material.  In  addition,  it  is  necessary  to  characterize  the  thermal  properties  of 
the  absorber  and  to  understand  how  the  heat  deposited  by  the  microwaves 
redistributes  within  the  absorber.  With  proper  understanding  of  these  three 
issues,  it  will  be  possible  to  associate  a  microwave  field  strength  with  an 
observed  temperature  from  the  infrared  camera. 
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2  Electromagnetic  Theory 


This  section  describes  the  mechanisms  through  which  electromagnetic 
energy  is  deposited  in  a  material  in  the  form  of  heat.  The  heat  deposition  is 
controlled  by  the  electromagnetic  field  strength,  spatial  distribution, 
orientation  and  the  electromagnetic  properties  of  the  material.  The  heat  flow 
and  final  temperature  profile  are  controlled  by  the  thermal  properties  of  the 
material  and  the  ambient  temperature.  The  emphasis  of  this  discussion  will 
be  on  the  electromagnetic  considerations. 


2.1  General  Description  of  the  Problem 

The  geometry  under  consideration  is  shown  in  Fig.  2-1.  A  slab  of 
microwave  absorber  is  entirely  surrounded  by  styrofoam.  (The  section  of 
foam  on  the  right  side  of  the  absorber  has  been  removed  so  that  the 
conduction  effects  can  be  illustrated.)  The  foam  provides  a  thermal 
insulation,  but  has  very  little  interaction  with  the  microwave  field.  For  the 
purposes  of  this  study,  the  effects  of  the  styrofoam  on  the  microwave  field 
will  be  ignored.  The  microwaves  interact  with  the  absorber  and  deposit 
energy  in  the  form  of  heat.  This  heat  can  redistribute  throughout  the 
absorber  by  means  of  the  thermal  conductivity.  In  addition,  the  heat  can 
leave  the  absorber  through  conduction,  convection  and  radiation  from  any  of 
its  sides.  Because  the  styrofoam  next  to  the  sides  of  the  absorber  will  heat  to 
the  same  temperature  as  the  absorber,  radiation  from  all  sides  except  the  top 
can  be  considered  as  negligible.  Because  the  styrofoam  touches  every  side 
except  the  top,  convection  from  all  sides  except  the  top,  can  be  ignored.  The 
total  amount  of  heat  lost  due  to  convection  is  proportional  to  surface  area. 

The  thickness  of  the  slab  will  typically  be  small  so  as  to  resemble  a  thin  sheet. 
Therefore,  the  heat  lost  through  the  thin  edges  can  be  ignored,  leaving  only 
the  heat  conducted  through  the  bottom  of  the  sheet  into  the  styrofoam  base. 
Finally,  the  sheet  of  absorber  is  placed  horizontally  so  that  the  convection  will 
not  distort  the  heating.  Analysis  of  the  thermal  parts  of  this  problem  will  be 
left  for  future  research. 


2.2  Electromagnetic  Analysis 

First,  the  electromagnetic  aspects  of  the  problem  outlined  above  will  be 
discussed.  This  will  include  considerations  of  how  heat  is  formed  in  a 
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Convection  Radiation 


material  in  the  presence  of  an  electromagnetic  field.  A  brief  discussion  will  be 
given  to  provide  physical  understanding  for  the  electromagnetic  parameters 
of  a  material.  Attention  will  be  given  to  how  the  orientation  of  the  field  lO  a 
thin  absorbing  screen  affects  the  heating.  The  objective  of  this  section  is  to 
adequately  describe  the  significant  factors  which  determine  how  much 
electromagnetic  energy  is  deposited  in  a  material. 


2.2.1  Conservation  of  Power  in  a  General  Medium 


The  constitutive  equations  relating  field  quantities  and  material 
parameters  in  the  frequency  domain  are 

D  =  £E 
B  =  pH 

J  =  oE+Ja 

In  the  last  equation,  Ja  represents  an  applied  or  source  volumetric  current. 
The  material  parameters  are  allowed  to  be  complex: 

e=  e'-je" 

0=0' -jo" 

This  allows,  for  example,  D  to  have  both  magnitude  difference  and  phase  lag 
with  respect  to  E  at  a  given  frequency.  For  this  work,  the  material  parameters 
are  considered  to  be  frequency  domain  quantities  which  may  be  functions  of 
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frequency.  Although  it  is  unlikely  that  any  real  materials  require  all  three 
quantities  to  be  complex,  this  theory  will  cover  the  most  general  case.  Also, 
in  order  to  allow  for  generality,  a  magnetic  current  term  will  be  included  only 
as  a  source,  thus  the  volumetric  magnetic  current  density  is 

M  =  Mg 

Ampere’s  law  is 

V  X  H  =  jcoD+  J 

Substituting  for  D  and  |  from  above  yields 

V  X  H  =  jco£E+  oE+  Ja 

Putting  in  the  complex  expressions  for  the  material  parameters  gives 

V  X  H  =  j(o{e' -  je")E+  (o'  -  ja")E+  Ja 

Forming  the  complex  conjugate  of  both  sides  and  separating  the  real  and 
imaginary  parts  of  the  E  multiplier  gives  the  result 

V  X  H*  =  -jco(e'  +  je")!’'  +  (o'  +  ja")E  V  Ja* 

V  X  h’'  =  [(o'  +  coe")  +  j(o"  -  tte')|E  V  Ja*  (2-1 ) 

We  will  return  to  this  result  shortly. 


Faraday's  law  is 


V  X  E  =  — jo>B  -  M 

Substituting  for  B  and  M  from  above  yields 

VxE=-j(DpH-Ma 

Putting  in  the  complex  expressions  for  the  material  parameters  gives 

VxE  =  -j(i)(p'-jp")H-Ma 

Now,  form  the  vector  inner  product  of  both  sides  with  the  complex  conjugate 
of  the  magnetic  field  H. 

(Vxe)  H*  =  -((op"  +  jo>p')H-  H’^-Ma-  H* 

The  left  side  can  be  simplified  by  means  of  the  vector  identity 

(Vx  A)  B=V-  (Ax  b)-!- A  (Vxb) 


which  yields 


V- 


ExH* 


+  E  Vx 


H’'j=-(cop"-(-jcop')H  H’^-Ma-  h’ 


Substituting  equation  2-1  into  this  yields 


ExH^j-t-E  [(o'-ha)e'')-tj(o''-(oe')]EVEla*  =  ^o)p''-(-jcop')H-  A’^-Ma-  iV 


The  identity  for  the  modulus  of  complex  numbers 


-  _  ♦  j-|2  2 

A  A  =|a|  =A^ 

can  be  used  to  simplify  some  of  the  dot  products.  Using  this  identity  and 
rearranging  gives 
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V-  ExH*j=-[(a'+coe'j  +  j(0''-(oe')]E^-E■%*-(co^"  +  )<o^')^^^-M  H* 
Now,  taking  the  volume  integral  of  both  sides  yields 

JJJV-  fExH“'ldv  =  -|JJ(o'  +  (oe'')E2+cog''H2dv-JJjMa-  H%Ela^clv- 

XT  \7 


jUJ(a''-a)e')E^+cog'H^dv 

V 

The  left  hand  side  can  be  reduced  by  the  divergence  theorem.  Rearranging 
yields 

-i|JjMa-  H’'+EVdv  =  “ffExH*j  ds  + 

^j|J(a'4-(oe")E2+o)g''H2dv+  ^2-2 


V 

" 

_v  V  ^  ^ 


Each  term  in  equation  2-2  can  be  identified.  This  yields  the 
conservation  of  power  equation: 

Ps  =  Pe+Pd+j2qWm-We) 

The  complex  power  supplied  by  the  sources  is 

Ps  =  -i|/jMaH’+EVdv 

V 

The  complex  power  exiting  the  volume  is 

Iff-  -♦V 
P  =_,£  ExH  ds 

The  real  power  being  dissipated  in  the  volume  is 

=  i  J  J  j(o'  +  Q)e")E^  +  cop"  d  V 

^  V 

The  time  averaged  stored  magnetic  energy  in  the  volume  is 

Wn-.=JJ|VH2dv 

V  ^ 

And  finally,  the  time  averaged  stored  electric  energy  in  the  volume  is 
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A  subtle  point  about  equation  2-3  should  be  emphasized;  it  is  necessary 
to  calculate  the  total  E  and  the  total  H  in  each  region  and  use  them  in  the 
equation.  In  the  case  of  plane  waves,  it  is  tempting  to  use  the  power  densities 
in  the  incident,  reflected  and  transmitted  waves  to  calculate  the  power 
absorbed.  However,  this  leads  to  erroneous  results.  A  simple  example  will 
serve  to  illustrate  this.  Consider  a  plane  wave  normally  incident  on  a 
boundary  as  shown  in  Fig.  2-2.  The  wave  reflects  from  the  boundary  with 
some  reflection  coefficient  F.  The  incident  electric  field  can  be  expressed  as 


Ei  =  zEoO  ^ 


where  y=a+jP  is  the  propagation  constant  with  a  the  attenuation  constant  in 
Nepers  per  meter,  P  the  phase  constant  in  radians  per  meter  and  the  field 
strength.  Thus,  the  power  density  in  the  incident  plane  wave  is 

211* 

The  power  density  in  the  reflected  plane  wave  is 


ri=^EixHi  =y 


2  ay 


Pr=^ErXH/=:-f 


l~Eo|' 

* 

211 


,2ay 


Thus  it  is  tempting  to  say  that  the  total  power  density  on  the  incident  side  is 
just  the  sum  of  these  two: 


ii  =  y 


2ti 


-2ay  _|j^2^2ay 


This  is,  however,  wrong. 


(2-4) 


The  total  electric  field  on  the  incident  side  is 

Et  =  Ei  +  Er  =  zEofe"^  +  Ve^  ^ 

V  ) 

The  total  magnetic  field  on  the  incident  side  is 

T1  V  ) 

Thus  the  total  power  density  on  the  incident  side  is 

Ft  =  -jEtx  Ht*  =  yl^  e"^“y  -i-  2jlm|rej^y  1  -|lj 


(2-5) 


Comparison  of  equations  2-4  and  2-5  reveals  a  missing  term  in  equation  2-4. 
This  term  is  important  in  that  it  maintains  the  conservation  of  power  for  all 
the  fields  in  this  example.  (Note  that  this  missing  term  does  not  contribute  to 
the  time-averaged  power  flow  because  the  time  averaging  involves  only  the 
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real  part  of  the  right-hand  side  of  equation  2-5.)  Therefore,  in  Section  3, 
where  plane  waves  incident  on  a  planar  absorber  are  considered,  rather  than 
differencing  the  incident  and  scattered  fields  to  obtain  the  power  dissipated,  as 
has  been  done  previously,  Pd  will  be  calculated  through  the  integral  indicated 
in  equation  2-3. 


Fig.  2-2  A  plane  wave  normally  incident  on  a  boundary  between 

2  materials 


2.2,2  Physical  Interpretation  of  Complex  Material  Parameters 

As  stated  at  the  beginning  of  Section  2.2.1,  the  electromagnetic  'material 
parameters  are,  in  general,  complex  values.  The  physical  interpretatu  it  of 
these  complex  values  will  now  be  discussed. 


2.2.2. 1  Polarization  Loss 

The  electric  flux  density  D  is  related  to  the  electric  field  E  and  the 
polarization  Pg  as 

D  =  £oE-i-  Pg 

For  a  linear  material,  the  polarization  is  linearly  dependent  on  the  electric 
field  as 

Pe  =  eoXeE 

where  a  dimensionless  quantity,  is  the  electric  susceptibility.  Thus,  the 
electric  flux  density  is  related  to  the  electric  field  as 

D  =  eo(l-l-Xe)E=eE 

Thus,  the  permittivity  of  the  material  is 

e  =  e'-je"  =  eo(l+Xe) 

Because  Eq  is  a  real-valued  constant,  Xg  must  be  complex  if  e"  ^  0.  Therefore, 
the  polarization  is  dependent  on  the  electric  field  in  a  complex  manner 
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I 


A  complex  relationship  between  two  quantities  implies  a  magnitude 
and  phase  difference.  Therefore^  the  polarization  must  be  out  of  phase  with 
an  applied  electric  field  when  e"  0.  From  equation  2-3,  it  can  be  seen  that 
the  term  o'  +  coe"  is  a  loss  term.  King  refers  to  this  quantity  as  the  "effective 
conductivity"  [2].  Thus,  the  imaginary  part  of  permittivity  leads  to  heating  of 
a  material.  Summarizing,  when  the  permittivity  of  the  material  is  complex, 
the  polarization  is  out  of  phase  with  the  electric  field  and  heat  is  generated  in 
the  material. 

A  physical  understanding  of  how  the  imaginary  part  of  permittivity  is 
similar  to  conductivity  can  be  gained  by  considering  an  electric  field  applied  to 
a  polarizable  material.  Consider  an  infinite  sheet  of  a  polarizable  material  in 
the  presence  of  a  uniform,  time  varying,  single  frequency  electric  field,  as 
shown  in  Fig.  2-3.  At  some  moment  of  time,  the  electric  field  is  maximum 
and  the  charges  are  distributed  as  shown  in  Fig.  2-3a.  Allowing  for  some 
phase  lag,  at  some  time  later  the  polarization  has  reached  a  maximum  and 
the  charges  are  as  shown  in  Fig.  2-3b.  At  this  moment,  however,  the  electric 
field  has  changed  from  that  in  Fig.  2-3a.  After  half  a  period  of  time  has 
elapsed  from  Fig.  2-3a,  the  electric  field  is  maximum  in  the  opposite  direction. 
The  charges  have  begun  to  realign  with  the  electric  field,  but  lag  somewhat 
behind,  as  shown  in  Fig.  2-3c.  Again,  allowing  for  some  phase  lag,  at  some 
time  later,  the  polarization  has  reached  a  maximum  in  the  opposite  direction 
from  Fig.  2-3b  and  the  charges  are  as  shown  in  Fig.  2-3d.  The  net  effect  is  that 
charges  are  in  motion.  These  charges  in  motion  are  equivalent  to  a  current. 
Thus,  the  polarizable  nature  of  the  material  has  led  to  a  "displacement 
current".  This  current  will  track  in  phase  with  the  electric  field  in  a  material 
whose  permittivity  is  entirely  real,  however,  the  phase  lag  in  a  material  with 
complex  permittivity  indicates  that  loss  is  occurring. 


+  +  4*  +  +  4*  +  4"  4"  4"  4"  4"  4- 4*  4"  4"  4- 4- 4" 


(a) 


(b) 


4-  4  4"  4"  4"  4  4"  4  4  4- 4"  4“  4"  4  4 


(c) 


(d) 


Fig.  2-3  Polarization  of  a  dielectric  sheet 


Magnetization  Loss 

The  magnetic  flux  density  B  is  related  to  the  magnetic  field  H  and  the 
magnetization  Pm  as 

B  =  fto(H+ fm) 

For  a  linear  material,  the  magnetization  is  linearly  dependent  on  the 
magnetic  field  as 

^m  “  Xm^ 

where  Xm'  ®  dimensionless  quantity,  is  the  magnetic  susceptibility.  Thus,  the 
magnetic  flux  density  is  related  to  the  magnetic  field  as 

B  =  go(l+Xm)H  =  gH 
Thus,  the  permeability  of  the  material  is 

p  =  p'-jp"  =  Po(l+Xm) 

Because  Po  is  a  real-valued  constant,  Xn\  niust  be  complex  if  p"  ^  0. 

Therefore,  the  magnetization  is  dependent  on  the  magnetic  field  in  a  complex 
manner. 

From  equation  2-3,  it  can  be  seen  that  the  term  p"  is  also  a  loss  term. 
Thus,  we  have  come  to  the  same  point  as  in  the  discussion  of  the  electric 
field.  Although  an  illustrative  example  showing  the  physical  explanation  of 
loss  due  to  complex  permeability  would  be  somewhat  more  complex  than 
that  for  permittivity,  the  conclusion  would  be  the  same,  namely,  the 
imaginary  part  of  permeability  results  in  a  phase  lag  between  the  magnetic 
field  and  the  magnetization  which  leads  to  heating  of  the  material. 


2.2.2.3  Conduction  Loss 

The  conduction  current  in  a  material  is  related  to  the  electric  field  as 

J  =  oE 

The  conductivity  is  allowed  to  be  complex  for  generality,  thus 

a  =  a'-ja" 

This  indicates  that  the  conduction  current  can  have  a  different  phase  from 
the  electric  field  in  a  material  where  o"  0.  This  is  not  the  case  for  the  type  of 
materials  considered  in  this  work.  However,  there  are  some  materials  which 
exhibit  complex  conductivity,  plasmas  for  example. 


The  materials  of  interest  in  this  work  will  typically  have  complex 
permittivities,  real  permeabilities  and  real  conductivities.  Equation  2-3  show.'^ 
that  the  two  parts  of  "effective  conductivity",  a  and  coc"  are  indistinguishable 
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as  far  as  their  effects  on  the  absorbed  power.  At  DC  only  a  affects  the  absorbed 
power.  At  high  frequencies  most  of  the  effective  conductivity  is  derived  from 
coe".  In  addition,  it  is  possible  to  use  an  effective  conductivity  which  groups 
all  the  electrical  loss  into  a  or  an  imaginary  permittivity  which  groups  all  the 
electrical  loss  into  coe".  Some  of  the  analyses  u.sed  in  this  work  use  cr  while 
others  use  coe";  the  form  chosen  is  that  which  seems  most  appropriate  to  the 
topic  of  discussion.  For  example,  the  measurement  of  material  parameters  at 
high  frequencies  outlined  in  Section  4  develops  equations  based  on  e".  Thus 
the  electrical  loss  at  high  frequencies  is  entirely  modeled  by  e".  Section  4 
develops  the  DC  material  measurement  on  the  basis  of  a.  To  add  a  to  the  wc" 
term  at  high  frequencies  is  not  allowable  because  the  coe"  term  already 
contains  all  the  lossy  characteristics. 


2.2.3  Power  Absorbed  Related  to  Field  Orientation 

From  the  previous  work,  the  power  dissipated  in  a  volume  due  to  the 
existing  electromagnetic  fields  in  that  volume  is 

Pd  =  -|  J /(o'  +  G)e")E^  +  cop"H^d  V 
^  V 

Consider  the  geometry  of  a  semi-infinite  slab  of  absorber  with  air  on  either 
side,  as  shown  in  Fig.  2-4.  It  is  preferable  to  talk  of  the  power  absorbed,  Pabs^ 
per  square  meter  (W/m^)  in  the  tangential  plane.  Thus, 

d 

Pabs  =  -j(('5'+<*")E2+fflti"H^dn  (2^6) 

0 

where  n  is  the  dimension  normal  to  the  boundary. 

Two  simplifying  assumptions  are  now  made.  First,  it  is  assumed  that 
the  absorber  thickness  d  is  very  small,  much  less  than  a  wavelength.  Second, 
it  is  assumed  that  the  material  parameters  are  such  that  the  absorption  is 
small.  As  a  result  of  these  assumptions,  the  fields  in  the  absorber  can  be 
approximated  as  constants  in  n. 

The  boundary  conditions  require  that  the  tangential  components  of  E 
and  H  be  continuous  at  both  boundaries.  From  the  thin,  low-loss 
assumptions  above,  it  is  possible  to  consider  all  tangential  E  components  to  be 
approximately  equal.  Thus, 

Etl  =  Et2=Et3  =  Pt 

where  the  subscript  refers  to  the  region  in  which  the  tangential  component 
exists.  The  same  applies  to  the  tangential  H  components.  Thus, 

^^tl  =  Ht2  =  Ht3  =  Ht 


K) 


Medium  1  Medium  2  Medium  3 


Fig.  2-4  Electromagnetic  fields  in  three  regions 

The  boundary  conditions  also  require  that  the  normal  components  of  D  and  B 
be  continuous  at  both  boundaries.  The  thin,  low-loss  assumptions  allow  the 
normal  D  components  in  all  regions  to  be  considered  approximately  equal. 
Thus, 

I^nl  ^n2  =  ^nS 

Expressed  in  terms  of  E  in  each  region  this  is 

eiEnl  =  e2En2  =  esEnS 

If  regions  1  and  3  are  both  vacuum  (or  air)  then  this  can  be  rewritten  as 

Enl  =  7^En2  =  En3  =  En 
^o 

Thus,  the  normal  component  of  the  electric  field  in  region  2  can  be 
approximated  as 

En2  =  En^ 

£2 

Similarly,  continuity  of  the  normal  component  of  B  yields 

Hn2  = 

^2 

Substituting  each  of  the  above  components  into  equation  2-6  and 
assuming  that  the  fields  do  not  vary  appreciably  in  n,  then  the  total  power 
absorbed  per  square  meter  in  region  2  is 


This  demonstrates  that  the  power  absorbed  from  the  normal  E  and  the 
normal  H  field  components  is  reduced  by  the  square  of  the  relative 
permittivity  and  permeability  of  the  absorber  material. 


Consider  an  absorbing  screen  which  is  designed  to  heat  from  the 


electric  field.  Most  electric  field  absorbers  are  likely  to  have  £r  2  3  or  4.  Thus, 
the  heating  due  to  the  normal  component  of  E  will  be  at  least  an  order  of 
magnitude  less  than  that  due  to  the  tangential  components. 

Therefore,  it  is  apparent  that  the  heating  of  a  thin,  low-loss,  absorbing 
screen  will  be  predominantly  from  tangential  field  components. 

The  amount  of  heating  which  arises  from  a  plane  wave  normally 
incident  on  a  semi-infinite  absorbing  sheet  will  be  considered  in  Section  3. 
This  particular  form  of  field  structure  will  be  used  to  design  an  absorber  with 
optimum  properties. 

3  Design  of  an  Absorbing  Screen 

The  typical  configuration  used  for  measuring  electric  fields  by 
thermographic  mapping  includes  a  thin  sheet  of  absorbing  material.  The 
design  objectives  for  this  material  are  that  it  exhibit  a  significant  heat  rise  for  a 
small  amount  of  incident  electromagnetic  energy  and  that  it  not  disturb  the 
fields  which  it  is  measuring.  These  two  objectives  are,  in  general,  in 
opposition.  Design  of  an  absorber  involves  a  trade-off  between  these  two 
desirable  properties. 

In  order  to  adequately  understand  the  trade-offs  in  the  design  of  an 
absorber,  it  is  desirable  to  be  able  to  predict  the  response  of  the  absorber  in 
some  situation  simple  enough  to  allow  for  a  closed-form  analysis.  One  such 
configuration  is  that  of  a  linearly  polarized  plane  wave  incident  on  a  semi¬ 
infinite  sheet  of  absorber.  This  configuration  will  be  analyzed  for  a  general 
case  and  then  reduced  to  the  case  of  normal  incidence.  The  results  from  the 
normal  incidence  analysis  will  be  used  to  evaluate  the  effects  of  the  absorber 
properties  on  the  heat  generated  in  the  absorber  and  the  amount  of  field 
disturbance. 


3.1  Plane  Waves  Incident  on  an  Absorbing  Screen 

The  geometry  to  be  considered  is  shown  in  Fig.  3-1.  The  incident  plane 
wave  is  linearly  polarized.  In  order  to  allow  for  generality  in  the  model,  both 
parallel  and  perpendicular  polarizations  (with  respect  to  the  plane  of 
incidence)  will  be  examined.  These  two  cases  yield  the  same  results  if  the 
angle  of  incidence  is  reduced  to  0  (normal  incidence).  The  incident  wave 


consisting  of  and  reflects  at  the  y=0  boundary  giving  rise  to  E^  and  ; 


it  also  transmits  through  the  boundary  giving  rise  to 


E'*'  and  H"^.  This  plane 
2  2 
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wave  reflects  at  the  y=d  boundary  giving  rise  to  and  it  also  transmits 

through  the  boundary  as  E  and  H  .  It  is  assumed  that  the  incident  field 

3  3 

strength  and  angle  of  incidence  are  known.  The  solution  for  the  unknown 
fields  and  angles  will  be  found  in  the  frequency  domain. 


3.1.1  Parallel  Polarization 

The  geometry  for  parallel  polarization  is  shown  in  Fig.  3-2.  The 
characteristic  impedance  and  propagation  constant  of  a  region  are  given  as 


Ya+ jcoe 

Y  =  jco^  |l+ ~  =  a -h  jP 
^  jcoe 


The  general  forms  for  the  electric  and  magnetic  fields  are 


Fig.  3-1  A  general  plane  wave  in  a  planar  3-region  geometry 
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X 

Fig.  3-2  A  parallel  polarized  plane  wave  in 

geometry 


a  planar  3-region 


1  ni 

E-  =  E^iysinei  +  zcose,)e-)P'(->'  ) 

^-^_fl.^-jP](-ycose,+zsine,) 

1  T|1 

E+  =  E+(-ysine2+ 

H+  =  ^  Je”''^2(y™592+zsine2) 

2  ^2 

E-  =  E-(ysine2  +  zcose2)e-^^<'>'  “*92+^sine2) 

H"  =  -^ie"^2(-y  cosSj+zsinez) 

2  T12 

E  =E  (-ysin93+zcose3)e-)'^3((>'^)“^®3+«in63) 

3  3 
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3  t13 

In  these  forms  of  the  fields,  regions  1  and  3  are  assumed  to  be  lossless. 
Furthermore,  it  is  assumed  that  regions  1  and  3  are  air  or  vacuum,  thus,  the 
intrinsic  impedances  and  propagation  constants  in  those  two  regions  are 

denoted  by  rio  and  jPo- 


Snell's  law  fixes  the  relationship  between  the  angles  as 

sin0i  Y2 


and  also 


sine2  jPo 

01=03 


Continuity  of  the  field  components  parallel  to  the  interfaces  yields  four 
boundary  conditions; 


+  E  =  E"^  +  E 
Iz  Iz  22  2z 


+  H  =  +  H 

lx  lx  2x  2x 


E  +E“  =E 
2z  2z  3z 


+H  =H 
2x  2x  3x 


ly=0,z=0 

y  =0,2=0 
ly=d,z=0 
y=d,z=0 


These  four  conditions  are  adequate  to  determine  E  ,  E'*’,  E  and  E  in  terms 

^  12  2  3 


of  E"^: 


E-  =  E-fc3 

^  ^  i-r2p2 

E'^=E'^ 


2  1  |_p2p2 

2  1  ,_r2p2 
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where: 


3  1  -j_p2p2 

11200502  -  TloCOS0| 
ti2Cos82+  iloCOsBi 
p  =  e''^2dcos92 


3.1.2  Perpendicular  Polarization 

The  geometry  for  perpendicular  polarization  is  shown  in  Fig.  3-3. 


y=0 


X 

Fig.  3-3  A  perpendicularly  polarized  plane 

region  geometry 


wave  in  a  planar  3- 


The  general  forms  for  the  electric  and  magnetic  fields  in  this  case  are 
given  as 

1  1 

H+  =  ^  (ysine,  -  zcos0i)e-)^  (ycos9,+zsine, ) 

1  Til 
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xe 


-jpl(-y  cos0i4-zsin0|) 


H~  =  — (ysin0|+  zcos0i)e 

1  711 


-jPi(-ycos0i+zsin0|) 


E'*'  =  E'^xe" 


-Y2(ycos02+zsin92) 


H^=-2.{ysin02-zcos02)e"^2(ycos02+zsin02) 

2  T12 

E"  =  E-5(e~^2(-y  cos02+zsin02) 

2  2 


H-  =  ^(ysine2+ 

2  T12 

E  =E  -e-jM(y-^)cos03+zsin03) 

3  3 

As  before,  Snell's  law  fixes  the  relationship  between  the  angles  as 

sin0i  Y2 
sin02  jPo 

and  also 

01-83 


Again,  continuity  of  the  field  components  parallel  to  the  interfaces 
yields  four  boundary  conditions: 


E"^  +  E  =  E"^  +  E 
lx  lx  2x  2x 


H'*’  +  H  =  H'*’  +  H 
Iz  Iz  2z  2z 


E+  +E^  =E^ 
2x  2x  3x 

+H~  =H 
2z  2z  3z 


ly=0,z=0 

Jy=0,z=0 
Jy=d,z=0 
y=d,z=0 
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3.2  Power  Absorbed  in  the  Screen 

Due  to  the  interaction  of  the  plane  wave  and  the  sheet  of  lossy 
material,  heat  will  be  deposited  in  the  material  as  described  in  Section  2.  The 
total  power  in  watts  absorbed  in  the  material  was  shown  to  be  related  to  the 
electromagnetic  field  by  the  relation 
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=  i|  j  j  (o'  +  coe")E^  +  cD^i"  d  v 

^  V 

For  plane  waves  it  is  possible  to  perform  the  indicated  integral  to  obtain  a 
closed-form  solution  for  P,^.  Because  the  intermediate  results  may  be  useful, 
they  will  be  given  here  as  steps  leading  to  the  closed-form  The  plane 
wave  is  assumed  to  be  normally  incident. 


The  total  electric  field  in  the  absorbing  region  is 

E  =)(E'^e~^2y+E-e+Y2y 
2  L  2  2 

The  magnitude  squared  of  the  electric  field  is 

2  2  2 

E  =  E^  =  E"^  e“^®2y  +  eH  e^“2y  +  2Rej E“^ E-^'^^yl 
2  2  2  2  1  2  2  J 

Integration  of  this  yields 


The  total  magnetic  field  in  the  absorbing  region  is 


H  =2-^e-'^2y+^e+Y2y 

2  \  r]2  112 

The  magnitude  squared  of  the  magnetic  field  is 


hJ  =Hi  = 


In  e 


2“2y  +  e;  2Re(E+Er*e-i2P2yl 

2  12  2  j 


Integration  of  this  yields 
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The  total  power  absorbed  per  square  meter  of  the  absorbing  material  is  thus 


The  baseline  design  for  the  screen  has  a  relative  permeability  ^ 

relative  permittivity  er=3,  a  conductivity  0=1  1/Qm,  a  thickness  d=0.003  m 
and  an  operating  frequency  of  f=3  GHz.  The  two  absorber  parameters  which 
were  most  easily  controllable  were  the  thickness  and  the  effective 
conductivity.  The  relative  permeability  and  permittivity  were  less 
controllable. 

In  the  following  series  of  plots  two  of  these  parameters  will  be  varied, 
while  all  others  will  stay  at  their  baseline  values.  The  plots  were  generated  by 
the  Mathematica  listing  shown  in  Appendix  A. 2.  In  these  plots,  the 
conductivity  was  used  rather  than  the  imaginary  part  of  the  relative 
permittivity  since  either  of  these  parameters  adequately  model  the  loss  term, 
as  discussed  in  Section  2.  An  alternative  method  for  calculations  based  on 
complex  permeability  and  complex  permittivity  is  given  in  the  Mathematica 
listing  in  Appendix  A.3. 

In  order  to  remove  the  effects  of  the  power  density  of  the  incident 
plane  wave  on  the  calculations,  the  power  absorbed  will  be  expressed  as  a 
percentage  of  the  incident  power  density,  i.e. 
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%Pabs  =  ^~-^'J00‘7« 

Mnc 

Fig.  3-4  shows  the  variation  of  %Pabs  varying  thickness  and 
conductivity  of  the  screen.  For  a  screen  with  a  thickness  of  4  cm  and  a 

conductivity  of  1  (Qm)'^  about  73%  of  the  incoming  planewave  power  can  be 
transformed  into  heat  in  the  screen.  It  should  be  noted,  however,  that  as  the 
material  is  made  thicker,  the  surface  temperature  is  reduced  due  to  thermal 
conduction  from  the  hot  surface  where  the  fields  are  more  intense  to  the  cool 
regions  inside  the  material  which  have  less  intense  fields.  Because  it  is  the 
surface  temperature  which  is  ultimately  of  interest,  it  is  not  practical  to  use 
thick  materials. 

Fig.  3-5  shows  the  variation  of  %Pabs  varying  frequency  and 
thickness.  This  Fig.  shows  that  the  absorber  can  be  constructed  to  be 
reasonably  insensitive  to  variations  in  frequency,  a  desirable  feature.  Fig.  3-6 
shows  the  variation  of  %Pabs  varying  relative  permittivity  and 
thickness.  This  shows  a  general  insensitivity  of  power  absorption  with 
relative  permittivity.  Fig.  3-7  shows  the  variation  of  %Pabs  varying 
relative  permittivity  and  frequency.  Fig.  3-8  shows  the  variation  of  %Pabs 
with  varying  conductivity  and  frequency.  Figs.  3-9  and  3-10  show  the 
variation  of  %Pabs  with  varying  thickness.  The  difference  between  the  two 
figures  is  in  the  range  of  thicknesses  considered. 
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Fig.  3-9  %Pabs  thickness 


Fig.  3-10  %Pabs  thickness  for  a  thin  absorber 


3.3  Power  Reflected  by  the  Screen 


Due  to  the  mismatch  between  the  intrinsic  impedance  of  the  air  and 
the  absorber  material,  some  power  will  be  reflected.  The  total  reflected  power 
in  watts/m2  is  given  as 


As  with  the  power  absorbed,  power  reflected  will  be  expressed  as  a  percentage 
of  the  incident  power  density,  i.e. 

%Prefl  =  X 100% 

Pine 

Fig.  3-11  shows  the  variation  of  %Prefl  with  varying  thickness  and 
conductivity.  As  the  conductivity  increases,  the  intrinsic  impedance  of  the 
material  is  less  well  matched  to  free  space  and  more  of  the  incident  wave  is 
reflected.  For  an  absorber  with  a  thickness  of  4  cm  and  a  conductivity  of  1 

(Dm)'^,  about  27%  of  the  power  is  reflected.  Together  with  the  73%  power 
absorbed,  this  accounts  for  the  power  flow  in  the  system. 

Fig.  3-12  shows  the  variation  of  %Prefl  with  varying  relative 
permittivity  and  thickness.  Figs.  3-13  and  3-14  show  the  variation  of  %Prefl 
with  varying  thickness. 
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Fig.  3-12  %Prefl  vs  relative  permittivity  and  frequency 
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Fig.  3-14  %Prefl  vs  thickness  for  a  thin  absorber 


3.4  The  Ratio  of  Absorbed  and  Reflected  Power 

It  is  generally  desirable  that  the  screen  absorb  much  more  power  than  it 


reflects.  The  effects  of  thickness  and  conductivity  on  the  ratio  of 
^Pabs/'^°Prefl  shown  in  Fig.  3-15.  The  peak  is  due  to  a  local  minimum  in 
Prefl-  Figs-  3-16  and  3-17  show  the  variation  of  %Pabs/^Fi-efi  with  varying 
thickness.  The  only  difference  between  the  two  figures  is  in  the  range  of 
thicknesses  considered. 


Fig.  3-15  %Pabs/^Fi-efl  vs  thickness  and  conductivity 
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Fig.  3-16  %Pabs/ '^F’refl  thickness 
%  Pabs/%  Prefl 


Fig.  3-17  %Pabs/‘^o^refl  vs  thickness  for  a  thin  absorber 


3.5  Absorber  Materials 


The  absorber  materials  chosen  for  fabrication  were  carbon-  and  ferrite- 
loaded  epoxies.  A  liquid  epoxy  system  is  particularly  useful  because  of  its 
ability  to  be  poured  into  molds  and  then  machined,  as  well  as  its  durability 
when  cured.  The  epoxy  resin  chosen  for  this  work  was  KPON-828  from  Shell 
Chemicals  Company.  A  variant  of  this  material  is  EPON-815  which  is  a 
thinned  version  of  the  828  resin.  If  further  thinning  of  the  resin  w'as 
necessary  in  order  to  prevent  excessive  thickening  at  high  loadings.  Shell 
BGE,  butyl  glyddyl  ether  was  used.  In  some  early  mixtures,  a  thin  epoxy 
called  Bisphenol-F  was  used.  There  are  no  advantages  of  this  over  a  thinned 
version  of  828  and,  thus,  its  use  was  discontinued.  The  curing  agent  for  all 
the  mixes  was  TETA  from  Pacific  Anchor.  All  of  these  materials  can  be 
obtained  from  EV  Roberts  &  Assoc,  in  Culver  City,  CA.  Cleanup  can  be 
accomplished  with  methyl  ethyl  ketone. 

The  carbon  loadings  used  Ketjen  EC  600  JD  from  Akzo  Chemicals  in 
Chicago,  IL.  This  is  a  very  low  density,  but  highly  conductive  carbon  black 
powder  which  resembles  the  toner  in  a  copy  machine.  The  ferrite  loadings 
used  Ferrite-50  from  Trans-Tech,  Inc.  of  Adamstown,  MD.  In  contrast  to  the 
carbon,  this  is  a  very  dense  powder. 

The  difficulties  in  working  with  the  two  loadings  are  at  opposite  ends 
of  the  viscosity  spectrum.  Because  the  carbon  has  a  very  low  density  (in  g/cc), 
mixtures  containing  more  than  about  2%  (by  weight)  carbon  in  815  epoxy  are 
prohibitively  thick.  If  a  higher  percentage  of  loading  is  desired  it  is  necessary 
to  thin  the  resin  with  BGE.  The  ferrite,  however,  has  a  very  high  density  and 
is  hard  to  keep  suspended  in  epoxy  without  falling  to  the  bottom  of  the 
mixture.  Mixtures  of  more  than  40%  ferrite  in  828  are  verv  difficult  to 
manufacture  with  good  homogeneity.  In  this  case,  it  would  be  necessary  to 
start  with  a  thicker  epoxy  resin.  Because  the  emphasis  of  this  work  was  on 
carbon-loaded  epoxies,  more  effort  was  devoted  to  thinning  to  obtain  good 
carbon  loadings  than  thickening  to  obtain  good  ferrite  loadings. 

In  order  to  prevent  air  bubbles  from  forming  in  the  epoxy  during 
mixing,  a  vacuum  mixer  was  used.  This  consisted  of  a  commercial  grade  food 
blender  with  a  stainless  steel  pitcher  and  a  reinforced  rubber  lid.  The  lid  was 
drilled  and  a  vacuum  fitting  installed.  The  blender  was  chosen  because  the 
mixing  blade  and  the  high  rpm's  aid  in  dispersing  the  loading  agent  in  the 
resin. 

3.6  Recipe  for  Preparation  of  Absorber 

The  epoxy  resins  must  be  used  in  correct  proportions  with  the 
hardener.  The  weight  ratio  of  TETA  to  828  is  142:1000,  while  the  ratio  oi 
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TETA  to  815  is  146:1000.  If  it  is  necessary  to  add  BGE  for  thinning,  then  extra 
TETA  should  be  added  in  the  ratio  208:1000,  TETA  to  BGE.  In  all  cases  the 
percentage  of  carbon  or  ferrite  loading  is  given  as  the  fraction  of  the  loading 
weight  to  the  final  mixture  weight. 

The  steps  taken  for  mixing  a  batch  of  loaded  epoxy  were; 

1)  Weigh  out  the  amount  of  loading  agent  and  put  aside 

2)  Weigh  out  the  amount  of  epoxy  resin 

3)  Add  thinner  to  the  resin  if  necessary 

4)  Add  half  of  the  hardener  to  the  results  of  3  and  mix  under 

vacuum  for  30  seconds 

5)  Add  one  quarter  of  the  loading  agent  and  mix  under  vacuum 

for  30  seconds 

6)  Repeat  step  5  twice  more 

7)  Add  the  second  half  of  the  hardener  to  the  results  of  6  and  mix 

under  vacuum  for  30  seconds 

8)  Add  the  final  quarter  of  the  loading  agent  and  mix  under 

vacuum  for  1  minute 

9)  Pour  into  a  mold  and  let  sit  at  room  temperature  for  at  least  5 

hours 

10)  Bake  the  epoxy  in  the  mold  for  2  hour  at  350°F. 

11)  When  cool,  remove  the  epoxy  from  the  mold 

Epoxies  hardened  with  TETA  are  best  hardened  when  the  mold  shape 
has  a  large  surface  area.  For  sample  weights  over  about  50  g,  the  curing  can 
become  exothermic  if  insufficient  surface  area  is  provided.  This  release  of 
heat  accelerates  the  cure  which,  in  turn,  releases  more  heat  with  the  result 
that  within  a  matter  of  a  few  minutes  the  material  will  harden.  This  means 
that  steps  4  to  8  must  be  accomplished  in  about  10  minutes. 

Glass  molds  were  tried,  but  proved  unusable  due  to  the  adhesion  of  the 
hardened  epoxy  to  the  glass.  Two  commonly  available  materials  which  do 
not  have  adhesion  problems  are  teflon  and  UHMW  (Ultrahigh  Molecular 
Weight  polyethylene).  Because  UHMW  was  far  less  expensive  than  teflon  it 
was  used  for  the  initial  molding  at  room  temperature.  The  hardened  resin 
was  removed  from  the  UHMW  after  5  hours  and  placed  between  teflon 
sheets  for  the  high  temperature  bake.  Both  of  these  materials  were  obtained 
at  Regal  Plastics  in  Colorado  Springs. 


35 


4  Measurement  of  the  Electromagnetic  Properties  of  the  Absorber 

One  area  which  has  been  lacking  in  previous  efforts  at  quantifying  the 
thermographic  mapping  of  electromagnetic  field  is  that  of  characterizing  the 
electric  and  magnetic  properties  of  the  absorbing  screen  at  microwave 
frequencies.  Some  of  the  previous  researchers  [3,4]  have  used  DC 
conductivity  measurements  to  characterize  the  loss  characteristics  of  the 
material.  Previous  attempts  at  measurement  of  the  permittivity  [3]  have 
suffered  from  unacceptably  large  variations  in  the  measured  results.  In  order 
to  enhance  the  understanding  of  the  material  parameters,  a  general  method 
will  be  developed  which  permits  determination  of  the  complex  permittivity 
and  permeability  at  microwave  frequencies  using  a  vector  network  analyzer. 
This  methodology  will  be  used  to  measure  well-known  materials  and  then 
those  fabricated  as  part  of  this  work.  Some  of  the  practical  considerations 
leading  to  more  av.curate  measurements  will  be  discussed. 

The  general  method  discussed  here  is  that  advocated  by  Weir  [5].  This 
method  has  been  shown  [6]  to  yield  accurate  results  when  the  loss  tangent  of 
the  material  is  greater  than  about  0.1.  The  method  of  measurement  is 
applicable  to  either  waveguide  or  coax,  but  the  derivation  here  will  focus  on 
coax.  First  the  voltage  and  current  in  a  coax  will  be  found  when  a  section  of 
lossy  dielectric  exists  between  2  lossless  regions.  These  will  then  be  converted 
to  S  parameters  to  permit  comparison  with  measurements  from  a  network 
analyzer.  A  conversion  between  S  parameters  and  material  parameters  will 
then  be  developed. 


4.1  Analysis  of  a  Section  of  Lossy  Dielectric  in  a  Coaxial  Airline 

Fig.  4-1  shows  a  section  of  lossy  dielectric  in  a  coaxial  airline.  The 

incident  voltage  and  current  wave  consisting  of  and  reflects  at  the  z=0 


boundary  giving  rise  to  and  I^;  it  also  transmits  through  the  boundary 
giving  rise  to  and  .  This  plane  wave  reflects  at  the  z=d  boundary  giving 


rise  to  V  and  I  ;  it  also  transmits  through  the  boundary  as  and  I'*’. 
2  2  3  3 


This 


geometry  and  that  for  plane  waves  given  in  Section  3  are  very  similar  and, 
hence,  lead  to  similar  forms  of  solution.  As  in  Section  3,  the  solution  for  the 
unknown  voltages  and  currents  will  be  found  in  the  frequency  domain. 
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The  characteristic  impedance  and  propagation  constant  of  a  section  of 
coaxial  cable  are  given  as 


1  2  3 


V^,I+ 

1  1 

v+  1+ 

2 ' 2 

V+,I+ 

3  3 

V  ,I 

1  1 

^2^2 

z=0  z=d 


Fig.  4-1  Lossy  dielectric  section  in  a  coaxial  airline 

In  these  two  equations  p  and  e  are  complex  and  represent  all  the  losses  from 
the  dielectric  section  of  the  coax.  All  the  effective  conductivity  is  expressed  in 

the  term  coe".  It  is  necessary  to  model  all  the  conductivity  in  this  manner 
because  the  network  analyzer  measurements  will  provide  two  complex 
quantities,  Sn  and  S21.  From  these  only  two  complex  quantities  can  be 

determined.  If  the  model  for  the  dielectric  included  complex  p  and  e  as  well 
as  a  conductivity  o,  the  number  of  unknowns  would  exceed  the  number  of 
measured  quantities  and  a  unique  solution  for  the  material  parameters 
would  not  be  possible.  Wall  losses  have  been  disregarded  in  these  forms  for  Z 
and  y. 


It  should  be  noted  at  this  p>oint  that  the  e"  which  is  measured  using  this 
method  is  not  the  same  as  that  which  was  part  of  the  electromagnetic  theory 
in  Section  2.  Let  us  call  the  lossy  part  of  the  permittivity  from  the 

measurement  being  considered  in  this  section  the  measured  imaginary 
permittivity.  The  relationship  between  this  quantity  and  the  lossy  electrical 
parameters  of  Section  2  is 

coe'm  =  cf  +  toe" 

Thus,  the  measured  permittivity  includes  both  the  material's  DC 
conductivity  and  the  material's  imaginary  permittivity.  However,  for  ease  of 
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notation,  the  measured  imaginary  permittivity  will  continued  to  be  referred 
to  as  e". 


The  general  forms  for  the  voltage  and  current  in  each  section  of  the 
coax  are 

11  1 

1  ZiV  1  1 

V  =V+e-Y2Z+Y-e+^2Z 
2  2  2 

2  Z2V  2  2  J 

Y  =  Y  ^  0“iP3(z"d) 

3  3 

3  Z3  3 

In  these  forms,  sections  1  and  3  are  assumed  to  be  lossless.  Furthermore,  it  is 
assumed  that  sections  1  and  3  have  air  as  the  dielectric,  thus,  the  characteristic 
impedances  and  propagation  constants  in  those  two  regions  are  denoted  by  Zq 

and  jPo- 

Continuity  of  the  voltage  and  current  at  the  boundaries  of  the  section 
of  dielectric  yields  four  boundary  conditioi\s: 


These  four  conditions  are  adequate  to  determine  V  ,  V^,  V  and  V"^  in 

^  12  2  3 

terms  of  V"^: 

1 
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V 


2 


V  =V' 
2  1 


V.'^  =  v 


I1 

i-p^ 

r  + 

v  ) 

J 1 

_r2p2 

/■+_ 

r+1 

’  1 

_r2p2 

.-rp^(r+ 1) 

1 

_r2p2 

H 

i-r^l 

.-1-  ' 

1  ) 

where: 


r= 


1  I_r2p2 
Z2-Zo 


Z2+Z0 

P  =  e-T^2d 

Assuming  that  the  left  side  of  the  coax  in  Fig.  4-1  is  connected  to  port  1 
of  a  network  analyzer  and  the  right  side  to  port  2,  the  measured  S  parameters 
are  related  to  the  voltages  as 


Sn  = 


V-  n 
1  _ 


i-p^ 


v+  i-r2p2 

1 


(4-1) 


S21 


V 

_  3 


+  pj^i-r^ 


v+  i-r2p2 

1 


(4-2) 


4.2  Relating  Material  Parameters  to  Measured  S-Parameters 

Equations  4-1  and  4-2  give  the  measured  S  parameters  in  terms  of  the 

reflection  coefficient  P  and  propagation  term  P.  A  very  tedious  inversion  of 
these  two  equations  yields 

r=x±i/x^-i  (4-3) 

Sii  +  S2i-r 

l-(Sll+S2i)r  (4-4) 
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where 


X  = 


2Sii 


(4-5) 


Because  of  the  symmetry  in  equations  4-1  and  4-2,  the  inversion  can  also  yield 

F=v±^^-1 


r= 


Sn  +  S2i-P 
l-(Sll+S2i)P 


where 


S21^~Sll^+  1 

2821 

The  form  in  equations  4-3  to  4-5  will  be  used  for  the  remainder  of  this  work. 

From  the  definitions  of  P  and  P,  the  material  parameters  can  be  found 
as  functions  of  P  and  P. 


=  y 


The  material  parameters  can  thus  be  found  as 

Pr2  =  xy 

er2  =  ^ 

X 


(4-6) 

(4-7) 

(4-8) 

(4-9) 


4.3  Data  Processing  for  the  Determination  of  Material  Parameters 

Determining  the  complex-valued  material  parameters  from  the 
measured  S  parameters  can  be  accomplished  by  application  of  equations  4-1  to 
4-9.  The  computer  code  used  to  perform  the  data  reduction  is  given  in 
Appendix  B.  There  are,  however,  two  issues  which  need  to  be  addressed 
before  this  data  reduction  can  occur.  These  are  the  subjects  of  the  next  two 
sections. 


4.3.1  Resolving  Multiple  Solutions  by  Velocity  Considerations 

Equation  4-7  has  a  multiplicity  of  solutions  due  to  the  complex  natural 
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logarithm.  Expanding  this  equation  yields 

1  nfe^^^e^^^l^ - J  (ad  +  j(pd  +  27i:n)) 

V  >'  (od^[io£o 

where  n  is  an  integer.  The  multiple  values  of  n  each  lead  to  different 
combinations  of  p  and  e  which  could  yield  the  measured  S  parameters.  That 
is,  the  mapping  of  S  parameters  to  p  and  e  is  not  1-to-l.  In  order  to  resolve 
this  ambiguity,  some  other  quantity  must  determine  which  value  of  n  is 
correct.  A  quantity  which  can  be  derived  from  the  measured  data  it.  the  group 

delay  of  the  signal.  Because  each  combination  of  p  and  £  leads  to  different 
predicted  group  delays  through  the  sample,  a  comparison  of  the  measured 
and  predicted  group  delays  will  determine  the  correct  value  of  n. 

The  group  delay  can  be  determined  from  the  measured  data  as 

5  / 

X(r  = - (< 

®  3to' 

where  <P  is  the  radian  phase  angle  of  the  term  P  in  equation  4-4.  Once  a 
value  of  n  is  chosen,  the  complex  material  parameters  for  that  choice  of  n,  p^ 
and  e^/  can  be  calculated.  The  group  delay  can  then  be  calculated  for  this 
choice  of  n  as 

_  d 

where 

^n  £n“}^n 

By  comparing  the  experimentally  determined  group  delay,  Xg,  and  the  group 

delay  for  a  particular  choice  of  n,  Xg^,  the  correct  value  of  n  can  be 

determined.  This  correct  value  of  n  then  yields  the  correct  values  of  p  and  e. 

Occasionally,  due  to  noise  in  the  data,  an  incorrect  choice  of  p  and  e  is  chosen, 
as  can  be  seen,  for  example,  in  the  spikes  in  Fig.  4-17. 


4.3.2  Network  Analyzer  Reference  Plane  Adjustment 


Another  modification  to  the  measured  data  involves  moving  the 
reference  plane  of  the  network  analyzer  to  each  side  of  the  sample.  Most 
network  analyzers  have  this  feature  built-in.  However,  for  this  work  a 
manual  adjustment  was  made  to  the  reference  planes.  Fig.  4-2  illustrates  the 
movement  of  the  reference  planes. 


The  network  analyzer  makes  measurements  of  S]|'  and  S21'  at  the 
points  at  which  it  was  calibrated.  However,  equations  4-1  through  4-5  use  the 
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values  just  at  the  ends  of  the  sample;  these  are  denoted  by  Sn  and  S21.  The 
conversion  from  the  measured  values  to  those  which  would  have  been 
measured  at  the  ends  of  the  sample  are 


Sn  =  Sii'ej^o^^ 
S2,  =  S2,'e»>^) 


Fig.  4-2  Network  analyzer  reference  plane  movement 


where 

Po  =  Wf^oEo 

It  should  be  noted  that  these  reference  plane  extensions  do  not  take  into 
account  losses  in  the  coax,  however,  for  a  good  quality  coax  these  losses  are 
small. 


4.4  Measurement  Equipment  and  Procedures 

Some  of  the  early  measurements  were  made  on  an  HP  8510A  network 
analyzer,  while  the  rest  were  made  on  a  HP  8720A.  A  full  two  port  calibration 
was  performed  prior  to  each  measurement  session.  The  reference  planes  after 
this  calibration  were  the  ends  of  the  cables  on  the  ports  of  the  analyzer.  A 
fairly  slow  sweep  rate  of  10  seconds  was  used  along  with  4  trace  averaging. 

The  data  was  captured  over  HPIB  and  processed  using  the  code  in  Appendix 
B. 


Two  very  important  quantities  are  the  reference  plane  offsets  a  and  b  in 
Fig.  4-2.  It  is  extremely  important  that  the  position  of  the  sample  in  the  coax 
be  known  to  at  least  ten  degrees  of  phase  or  less.  As  an  example,  a  1  mm 
uncertainty  in  the  position  of  the  sample  at  6  GHz  is  7.2  degrees  of  phase,  but 
at  10  GHz  it  is  12  degrees.  This  is  critical  because  the  complex  material 
parameters  are  determined  from  complex  S  parameters.  A  variation  in  phase 
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of  the  S  parameters  can  cause  large  variations  in  the  material  parameters. 

This  partly  explains  the  instability  in  the  data  at  higher  frequencies.  Most 
measurements  with  the  network  analyzer  involve  magnitude  only  data.  The 
need  for  accurate  phase  data  makes  the  measurement  more  difficult.  A  jig  to 
precisely  position  the  sample  in  the  coax  is  quite  helpful  at  reducing  the 
uncertainty  in  the  sample's  position.  The  measurements  showed  good 
repeatability  when  the  alignment  jig  was  used. 

Because  the  measurement  is  being  made  at  the  end  of  two  to  three  foot 
long  teflon  dielectric  cables,  a  short  and  low  loss  sample  may  not  cause  a 
significant  change  in  the  S  parameters.  Longer  samples  have  more  effect  on 
the  S  parameters  and,  thus,  provide  better  data.  The  same  is  true  of  higher 
loss  samples.  Hewlett  Packard  believes  that  this  method  is  accurate  if  the  loss 
tangent  is  greater  than  about  0.1  [6].  Although  no  attempt  was  made  to  verify 
this  claim,  the  results  do  seem  better  for  the  higher  loss  samples. 

One  final  note  about  the  measurements  is  in  order.  Due  to  the 
uncertainty  in  the  measurement,  the  data  will  typically  yield  bad  values 
when  the  sample  is  a  multiple  of  a  quarter  of  a  wavelength.  (This  is  a  quarter 
of  a  wavelength  in  the  sample  material,  not  free  space.)  The  remedy  for  this 
is  to  make  the  sample  short  enough  that  the  first  quarter-wavelength 
frequency  is  above  the  measurement  band.  Notice  that  this  is  in  opposition 
to  the  desire  for  a  long  sample  size  stated  above.  In  order  to  resolve  these 
opposing  constraints,  the  sample  length  for  the  loaded  epoxies  was  chosen  to 
be  as  long  as  possible  while  still  staying  just  below  the  frequency  of  the  first 
discontinuity.  This  was  not  done  for  the  air  or  nylon  samples  and,  therefore, 
plots  for  these  materials  display  the  discontinuities. 


4.5  Material  Measurement  Results 


Figs.  4-3  and  4-4  show  the  complex  material  parameters  for  an  air  filled 
10  cm  coax.  The  data  are  within  2%  of  the  expected  values  except  at  the 
quarter  wavelength  frequencies.  The  imaginary  parts  should  be  zero,  but  are 
not.  This  could  be  due  to  experimental  uncertainty  or  to  incorporation  of 
wall  losses  as  dielectric  losses,  since  the  theory  does  not  account  for  wall 
losses.  As  is  normally  the  case  for  this  method,  the  lower  frequencies  are 
more  accurate  than  the  higher  frequencies. 
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Air 


Frequency 

Fig.  4-3  Real  permittivity  and  permeability  of  air 


Air 


Frequency 

Fig.  4-4  Imaginary  permittivity  and  permeability  of  air 
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4.5.1  Comparison  with  Other  Measurements 


Figs.  4-5  and  4-6  show  the  real  and  imaginary  parts  of  the  relative 
permittivity  and  permeability  for  teflon.  Von  Hippel  [7]  gives  the  value  of 

ei.=2.1-j3.15xl0"^  at  3  GHz.  Because  teflon  is  essentially  non-magnetic, 

As  discussed  above,  the  measurement  of  the  imaginary  parts  of  permittivity 
and  permeability  are  not  accurate  due  to  the  low-loss  nature  of  the  material. 
However,  the  real  parts  are  in  good  agreement  with  the  previously  published 
values. 

Figs.  4-7  and  4-8  show  the  real  and  imaginary  parts  of  the  relative 
permittivity  and  permeability  for  nylon.  Von  Hippel  gives  the  value  of 

er=3.03-j4.22xl0'3  at  3  GHz.  As  with  teflon,  nylon  is  essentially  non-magnetic, 

thus, 


4.5.2  Results  for  Loaded  Epoxies 

Each  of  the  plots  in  this  section  show  the  variation  of  the  real  and 
imaginary  permittivity  and  permeability  with  frequency.  In  addition,  an 
effective  conductivity  a^-cocQE"  has  been  plotted  with  the  imaginary  parts. 

Figs.  4-9  and  4-10  show  the  measured  material  parameters  for  a  20% 
ferrite  loading.  The  recipe  for  this  mixture  is  shown  in  Table  4-1. 
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O.Oe+00  l.Oe+09  2.0e+09  3.0e+09  4.0e+09  5.0e+09 

Frequency 

Fig.  4-7  Real  permittivity  and  permeability  of  nylon 


O.Oe+00  l.Oe-i-09  2.0e+09  3.0e+09  4.0e+09  5.0e+09 

Frequency 

Fig.  4-8  Imaginary  permittivity  and  permeability  of  nylon 


Ferrite  20% 


Frequency 

Fig.  4-10  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  20%  ferrite  loaded  epoxy 


Table  4-1  20%  ferrite  recipe 


Additive 

grams 

EPON-828 

50 

Bisphenol-F 

20 

Ferrite-50 

20.09 

TETA 

10.38 

Figs.  4-11  and  4-12  show  the  measured  material  parameters  for  a  30% 
ferrite  loading.  The  recipe  for  this  mixture  is  shown  in  Table  4-2. 


Table  4-2  30%  ferrite  recipe 


Additive 

grams 

EPON-828 

100 

Ferrite-50 

48.9 

TETA 

14.2 

Figs.  4-13  and  4-14  show  the  measured  material  parameters  for  a  40% 
ferrite  loading.  The  recipe  for  this  mixture  is  shown  in  Table  4-3. 
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Table  4-3  40%  ferrite  recipe 


Additive 

EPON-828 

200 

Ferrite-50 

97.9 

TETA 

28.4 

Figs.  4-15  and  4-16  show  the  measured  material  parameters  for  a  1.7% 
carbon  loading.  The  recipe  for  this  mixture  is  shown  in  Table  4-4. 


Ferrite  30% 


Frequency 

Fig.  4-11  Real  permittivity  and  permeability  of  30%  ferrite 

loaded  epoxy 


50 


0.30 


Ferrite  30% 


Frequency 

Fig.  4-12  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  30%  ferrite  loaded  epoxy 


t 

I 


51 


Ferrite  40% 


Frequency 


Fig.  4-14  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  40%  ferrite  loaded  epoxy 
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Carbon  1.7% 


O.Oe+00  l.Oe+09  2.0e+09  3.0e+09  4.0e+09  5.0e+09 

Frequency 

Fig.  4-16  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  1.7%  carbon  loaded  epoxy 


Table  4-4  1.7%  carbon  recipe 


Additive 

grams 

EPON-828 

24.125 

Bisphenol-F 

20.76 

EC  600  P 

0.875 

TETA 

5.71 

Figs.  4-17  and  4-18  show  the  measured  material  parameters  for  a  4% 
carbon  loading.  The  recipe  for  this  mixture  is  shown  in  Table  4-5. 


Table  4-5  4%  carbon  recipe 


Additive 

EBB 

EPON-815 

150 

BGE 

6 

EC600P 

7.5 

TETA 

23.6 

The  desired  parameters  for  the  absorbing  screen  were  o=l  (wc"=6),  f'=3 
and  p'=l  at  3  GHz.  The  4%  carbon  results  come  close  on  all  but  the  c'.  This  is 
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due  to  the  carbon  increasing  both  the  real  and  imaginary  parts  of  e. 

Figs.  4-15  and  4-16  show  the  measured  material  parameters  for  a  2.3% 
carbon  and  22.5%  ferrite  loading.  The  recipe  for  this  mixture  is  shown  in 
Table  4-6.  This  mixture  ratio  gives  the  percentage  of  loadings  disregarding  the 
weight  from  the  other  loading  component. 


Carbon  4% 


Frequency 

Fig.  4-17  Real  permittivity  and  permeability  of  4%  carbon  loaded 

epoxy 
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Carbon  4% 


Frequency 

Fig.  4-18  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  4%  carbon  loaded  epoxy 
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Carbon  2.3%  Ferrite  22.5% 


Frequency 


Fig.  4-20  Imaginary  permittivity,  permeability  and  effective 
conductivity  of  2.3%  carbon-22.5%  ferrite  loaded  epoxy 

Table  4-6  2.3%  carbon-22.5%  ferrite  recipe 


Additive 

grams 

EPON-815 

150 

EC  600  JD 

4 

Ferrite-50 

50 

TETA 

21.9 

Table  4-7  contains  a  summary  list  of  material  parameters  at  selected 
frequencies  for  the  samples  above. 
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Table  4-7  Material  parameters  for  prepared  samples 


Sample 

Frequency 

(GHz) 

^Bi 

20%  ferrite 

3.0 

i£iwgsm«ii 

30%  ferrite 

3.0 

■WiteglRgJtl 

40%  ferrite 

3.0 

1.7%  carbon 

1.0 

4%  carbon 

3.025 

■ilgrf;5j6jkfcfci 

2.3%  carbon- 
22.5%  ferrite 

3.0 

1.115-j.l425 

4.6  DC  Conductivity  Measurements 


The  measurement  of  conductivity  using  a  standard  method  is 
illustrated  in  Fig.  4-21.  The  total  resistance  of  the  block  of  material  is  given  as 


R  = 


The  conductivity  is,  thus,  given  as 


L 

CTwh 


L 

CT  = - 

Rwh 

The  carbon  4%  sample  measurement  had  a  DC  conductivity  of  0.066  (1/Qm). 
Fig.  4-22  shows  the  extrapolation  of  the  high  frequency  measurement  of  the 

effective  conductivity  (-toeoe")  down  to  a  DC  value  by  means  of  a  straight  line. 


This  predicts  a  DC  conductivity  of  about  0.05  1/Dm,  in  good  agreement  with 
the  measured  value.  Recall  that  the  measured  imaginary  permittivity 
includes  both  the  material's  DC  conductivity  and  the  material's  imaginary 
permittivity,  as  explained  in  Section  4.1.  Thus,  the  extrapolation  of  the 
measured  imaginary  permittivity  down  to  DC  should  yield  a  value  close  to 
the  measured  DC  conductivity,  as  is  the  case  in  Fig.  4-22. 
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Frequency 


Fig.  4-22  Extrapolated  effective  conductivity  for  4%  carbon 

loaded  epoxy 
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Appendix  A  Mathematica  Equations  for  Electromagnetic  Analyses 

This  section  lists  the  the  equations  used  in  Mathematica  for  the 
generation  of  the  plots  in  Section  3. 


A.l  Absorber  Equations 

The  following  input  can  be  used  to  calculate  the  power  flow  in  and 
around  a  semi-infinite  slab  of  absorber. 

eps0=8.854H0'^-12; 

mu0=4*PiH0^-7; 

etaair=N[Sqrt[muO/epsO]]; 

pincside[freq_,epsrdiel_,d_,sigma_]:= 

Block[ 

{omega,gammadiel,alpha,beta,GtadiGl,p,gamma, 

erefl,powerpersqm), 

omega=N[2*Pi*freq]; 
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gammadiel=N[I*omega*Sqrt[muO*epsO*epsrdiel]* 

SqrtIl+sigma/(romega*epsO*epsfdiel)JJ; 

alpha=Re[gammadiel]; 

beta=Im[gammadiel]; 

etadiel= 

N[Sqrt[muO/(epsO*epsrdiel+sigma/(romega))}j; 

p=N[E'^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(ctadiel+etaair)]; 

erefl=N[gamma*(l-p'^2)/(l-(gamma'^2)(p^2))J; 

powerpersqm=l  /(l^etaair)* 

(l-Abs[erefl]'^2+2*PIm[erefl]); 

Return[N[powerpersqm]]; 


pabs[freq_,epsrdiel_,d^sigma_]:= 

Block[ 

{omega,gainmadiel,alpha,beta,etadiel,p,gamma, 

e2plus,e2minus,powerpersqm}, 

omega=N  [2’^Pi*f  req] ; 

gammadiel=N[I*omega*Sqrt[muO*epsO*epsrdiel]* 

Sqrt[l+sigma/(I*omega*epsO*epsrdiel)]]; 

alpha=Re[gammadiel]; 

beta=Im[gammadiel]; 

etadiel= 

N[Sqrt[muO/(epsO*€psrdiel+sigma/(I*omega))]]; 

p=N[E^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(etadiel+etaair)]; 

e2plus=N[(gamma+l)/(l-(ganuna'^2)(p''2))]; 
e2minus=N  [(-gamma*p''2)*(gamiT\  a+ 1 )  / 
(l-(gamma^2)(p'^2))]; 

powerpersqm=sigma/4* 

((Abs[e2minus]^2)/aIpha*(E^(2*alpha*d)-l)- 
( Abs  [e2plus]^2)  /  alpha*(E^  (-2*alpha*d)- 1 ) + 
(2/beta)*Re[e2pIus*Conjugate[e2minus]* 

P(EA(.2*l’^beta*d)-l)]); 


Return[N  [powerpersqm]]; 

1; 


ptrans(freq_,epsrdiei_,d_,sigma_]:= 


Block[ 

{omega,gammadiel,alpha,beta,etadiel,p,gamma, 

etrans,powerpersqm}, 

omega=N[2*Pi*freq]; 

gammadiel  =N  [I*omega*Sqr  t[muO*epsO*epsrdiel  ]  * 
Sqrt[  1  +sigma  /  (I*omega*epsO*epsrdiel)  ]  ]; 
alpha=Re[gammadiel]; 
beta=Im[gammadiel]; 
etadiel= 

N[Sqrt[muO/(epsO*epsrdiel+sigma/(I*omega))]]; 

p=N[E^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(etadiel+etaair)]; 

etrans=N[p*(l-gamma''2)/(l-(gamma'^2)(p^2))]; 
powerpersqm=(  Abs[etTans]^2)  /  (2*etaair); 

Return[N[powerpersqm]]; 


eenergystored  [freq_,epsrdiel_,d_,sigma_] : = 

Block[ 

{omega,gammadiel,alpha,beta,etadiel,p, gamma, 
e2plus,e2minus,eener  gy  persqm} , 

omega=Nl2*Pi’*^freq]; 

gammadiel=N[Pomega*Sqrt[muO*epsO*epsrdiel]* 
Sqrt[l  +sigma/  (Pomega*epsO*epsrdiel)  ]  ]; 
alpha=Re[gammadiel]; 
beta=Im[gammadiel]; 
etadiel= 

N[Sqrt[muO/(epsO*epsrdiel+sigma/(I*omega))]]; 

p=N[E'^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(etadiel+etaair)]; 

e2plus=N[(gamma+l)/(l-(gamma^2)(p'^2))]; 

e2minus=N[(-gamma*p''2)’^(gamma+l)/ 

(l-(gamma^2)(p'^2))]; 

eenergypersqm= 

I*omega*epsrdiel’^eps0/4* 

((Abs[e2minus]^2)/alpha*(E^(2*alphaM)-l)- 

(Abs[e2plus]^2)/aIpha’^(E''(-2*aIpha*d)-l)+ 

(2/beta)*Re[e2plus’^Conjugate[e2minus]* 

I*(E^(-2*rbetaM)-l)]); 


Return  [N  [eenergy  pers  qm  ]  ]; 


henergystored[freq_,epsrdiel_,d_,sigma_]:= 

Block[ 

{omega,gammadiel,alpha,beta,etadiel,p,gamma, 

e2plus,e2minus,henergypersqm}, 

omega=N[2*Pi*freq]; 

gammadiel=N[I*omega*Sqrt[muO*epsO*epsrdiel]* 

Sqrt[l+sigma/(I*omega*epsO*epsrdiel)]l; 

alpha=Re[gammadiell; 

beta=Im[gammadiel]; 

etadiel= 

N[Sqrt[muO/ (epsO*epsrdiel+sigma/(I*omega))]]; 
p=N[E^(-gammadielM)l; 
gamma=N[(etadiel-etaair)  /  (etadiel+etaair)  ] ; 

e2plus=N[(gamma+l)/(l-(gan\ma''2)(p'^2))]; 

e2minus=N[(-gamma*p'^2)*(gamma+l)/ 

(l-(gamma'^2)(p'^2))]; 

henergypersqm= 

Pomega*muO/(4*Abs[etadiel]''2)* 

(( Abs[e2mmus]  ^^2)  /  alpha*(E^(2*alpha’‘d)-l)  - 

(Abs[e2plus]''2)/alpha*(E''(-2*alpha*d)-l)- 

(2/beta)’*Re[e2plus*Conjugate[e2minus]* 

r(E^(-2=»PbetaM)-l)]); 

Return[N[henergypersqm]]; 


A.2  Absorber  Design  Equations 

The  following  equations  were  used  to  generate  the  data  for  the  plots 
section  3. 

eps0=8.854n0''-12; 

mu0=4’»Pi*10''-7; 

pinc=1.0/377.0/2.0; 

gammaair=N[I*omega’^qrt[muO*epsO]]; 

etaair=N[Sqrt[muO/epsO]]; 


pabs[freq_,epsrdiel_,d_,sigma_]  := 


BIock[ 

{omega,gammadiel,alpha,beta,etadiel,p,gamma, 

e2plus,e2minus,powerpersqm), 

omega=N[2*Pi*freq]; 

gammadiel=N[I*omega’*Sqrt[muO*epsO*epsrdiel]* 
Sqrtll  +sigma  /(I*omega*epsO*epsrdiel)]]; 
alpha=Re[gammadiell; 
beta=Im[gammadiel]; 
etadiel= 

N[Sqrt[muO/(epsO*epsrdiel+sigiria/(I*omega))]]; 

p=N[E'^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(etadiel+etaair)]; 

e2plus=N[(gamma+l)/(l-(gamma'^2)(p'^2))]; 

e2minus=N[(“gamina*p''2)*(gamma+l)/ 

(l-(gamma‘^2)(p^2))]; 

powerf)ersqm=N  [sigma/  (4*alpha)* 

((Abs[e2minus]'^2)*(E‘^(2*alpha*d)-l)- 

(Abs[e2plus]^2)’*(E''(-2*alpha*d)-l))+ 

sigma/(2*beta)*Re[e2plus*Conjugate[e2minus]* 

I*(EA(.2*l*beta*d)-l)]]; 


Re  turn[N  [  Abs[  po  werpersqm]  ]  ]; 

]; 

erefl[freq_,epsrdieI_,d_jSigma_]:= 

Block[ 

{omega,gammadiel,alpha,beta,etadiel,p,gamma}, 

omega=N[2*Pi*freq]; 

gammadiel=N[I*omega*Sqrt[muO*epsO*epsrdiel]* 
Sqrt(l  +sigma/ (romega*epsO*epsrdiel)]]; 
alpha=:Re[gammadiel]; 
beta=Im[gammadiell; 
etadiel= 

N[Sqrt[muO/(epsO’^epsrdiel+sigma/(I*omega))]]; 

p=N[E'^(-gammadielM)]; 

gamma=N[(etadiel-etaair)/(etadiel+etaair)]; 

Return[N[Abs[gamma*(l-p'^2)/(l-(gamma'^2)(p'^2))]]]; 

]; 

prenpercGnt[freq_,epsrdiel_jd_,sigma_];= 
(GrGfl[freq,epsrdiel,d,sigma]  '^2)*1 00.0; 
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pabspercent[freq_,epsrdiel_,d_,sigma_]:= 

pabs[fTeq,epsrdiel,d^igma]  /pinc*l  00. 0; 


A.3  Absorber  Design  Equations  for  Complex  Material  Parameters 

The  following  equations  are  alternate  forms  to  those  in  section  A. 2. 
The  forms  here  allow  for  a  specification  of  complex  permittivity  and  complex 
permeability. 

eps0=8.854*10''-12; 
mu0=4»Pin0^-7; 
pinc=l  .0/377.0/2.0; 

gammaair=N[Pomega*Sqrt[mu0*eps0]]; 
etaair=N{Sqrt[muO /epsO]  ]; 

pabs[freq_^epsrdiel_,murdiel_,d_]:= 

Block  [ 

{mudiel,epsdiel,omega,gammadiebalpha,beta,etadiel, 

p,gamma,e2plus,e2minus,powerpersqm), 

omega=N[2*Pi*freq]; 

mudiel=N[murdiel*muO]; 

epsdiel=N[epsrdiel’^epsO]; 

gammadiel=N[I*omega’^Sqrt[mudiePepsdiel]]; 

alpha=Re[gammadiel]; 

beta=Im[gammadiel]; 

etadiel=N[Sqrt[mudiel/epsdiel]]; 

p=N{E'^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/{etadiel+etaair)]; 

e2plus=N[(gamma+l)/(l-(gamma'^2)(p''2))]; 

e2minus=N[(-gamma*p^2)*(gamma+l)/ 

(1 -(gamma^2)(p‘^2))]; 

powerpersqm=N[omega/4* 

((Abs[e2minus]^2)*Im[epsdiel]/alpha* 

(E^(2*alpha*d)-1)- 
( Abs[e2plus]  '^2)*Im[epsdiel]  /  alpha* 

(E^(-2*alpha*d)-l)+ 

(2*Im[epsdiel]/beta)* 

Re[e2pIus*Conjugate[e2minus]* 

I*(E'^(-2*rbeta*d)-l)]+ 

(Abs[e2minus]^2)*Im[mudicI]  / 
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(alpha*  Abs[etadiell'^2)* 

(E^(2*alpha*d)*l)- 
(Abs[e2plus]^2)*lmlmudiel]  / 
(alpha*Abs[etadiel]''2)* 

(E'^(-2*aIpha*d)-l)- 

(2*Im[mudiel]/(beta*AbsIetadiel]'^2))* 

Re[e2plus*Conjugate[e2minus]* 

I*(E'^(-2*I*beta*d)-l)])]; 

Retum[N[Abs[powerpersqm]]]; 

1- 

erefl  [  freq_,epsrdiel_,murdiel_,d_]  ;= 

Block[ 

{mudiel,epsdiel,omega,gammadiel,etadiel,p,gamma}, 

omega=N[2*Pi*freql; 

mudiel=N[murdiel*muO]; 

epsdiel=N  [epsrdiel*epsO]; 

gammadieI=N[I*omega*Sqrt[mudiel*epsdiel]]; 

etadiel=N[Sqrt[mudiel/epsdiel]]; 

p=N[E'^(-gammadiel*d)]; 

gamma=N[(etadiel-etaair)/(etadiGl+etaair)]; 

Return[N[Abs[gamma*(l-p'^2)/(l-(gamma'^2)(p^2))]]]; 

]; 

preflpercent[freq_,epsrdiel_,miirdiel_,d_]:= 

(erefl[freq,epsrdiel,murdiel,d]'^2)*100.0; 

pabspercent[freq_,epsrdiel_,murdiel_,dj:= 

pabs[freq,epsrdiel,murdiel,d]  /pinc*l  00.0; 
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Appendix  B  Computer  Processing  of  Material  Parameter  Data 

#define  AGREEMENT  0.2 

mainO 

{ 

register  int  loop,loop2; 

int  num_points,agreement,branch; 

char  buffer_arr[80]/buffer=(char*)&buffer_arr; 

double  d,portl_offset,port2_offsetdum_double,slope, 

exper_group_delay,branch__group_delay,dissipation_ratio, 
phase_sum,freq_sum,freq2_sum,product_sum; 
double  *freq/p_phase; 

complex  mu_r,epsilon_r,temp_sll,temp_s21,dum_complex; 
complex  *sll,*s21/p; 

FILE  *input_filesl  l/input_files21,*output_file; 


freq=NULL; 

p_phase=NULL; 

sn=NULL; 

s21=NULL; 

p=NULL; 

input_filesll=NULL; 

inputJiles21=NULL; 

output_file=NULL; 

freq=(double*)calloc((unsigned)405,(unsigned)sizeof(double)); 
if(freq==NULL)  { 

printfC'Error  allocating  freq  arrayXn"); 
goto  stop_program; 

} 

p_phase=(double’^)calloc((unsigned)405,(unsigned)sizeof(double)); 
if(p_phase==NULL)  { 

printfC'Error  allocating  p_phase  arrayXn"); 
goto  stop_program; 

} 

sll=(complex*)calloc((unsigned)405,(unsigned)sizeof(complex)); 
if(sll==NULL)  { 

printfC’Error  allocating  sll  array \n"); 
goto  stop_program; 

} 

s21=(complex*)calloc((unsigned)405,(unsigned)sizeof(complGx)); 
if(s21==NULL)  { 

printfC'Error  allocating  s21  arrayXn"); 
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goto  stop_program; 

} 

p=(complex*)calloc((unsigned)405,(unsig;ned)sizeof(complex)); 
if(p==:NULL)  { 

printfC’Error  allocating  p  arrayXn"); 
goto  stop_program; 

) 

printfC'Input  file  with  sll  data:  "); 
buffer=cgets(buffer); 
input_filesll=fopen(buffer/'r'’); 
if(input_filesn==NULL)  { 
printfC’File  not  found  \n"); 
goto  stop_program; 

} 

printfC'Input  file  with  s21  data:  "); 
buf  fer=cgets(buf  fer) ; 
input_files21=fopen(buffer/'r"); 
if(input_files21==NULL)  { 
printfC'File  not  foundXn"); 
goto  stop_program; 

} 

printfC'Output  file  :  "); 
buffer=cgets(buffer); 
output_file=fopen(buffer/'w"); 
if(output_file==NULL)  { 

printfC'File  not  found  \n"); 
goto  stop_program; 

} 

printfC'Thickness  of  Material  (cm):  "); 
buffer=cgets(buffer); 
sscanf(buffer/'%lE",&d); 
d=d/100.0; 

printfC'Distance  from  port  1  ref  plane  to  sample  (cm):  "); 
buffer=cgets(buffer); 
sscanf(buffer,"%lE",&portl_offset); 
portl_offset=portl_offset/100.0; 

printfC'Distance  from  port  2  ref  plane  to  imple  (cm):  "); 
buffer=cgets(buffer); 
sscanf(buffer,"%lE",«&port2_offset); 
port2_offset=port2_offset/ 100.0; 
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num_points=0; 

while(!feof(input_filesll)  &&  !feof(input_files21))  { 
if(fgets(buffer,80,input_filesl  1)!=NULL)  { 
i\um_points++; 

sscanf(buffer,"%lE%lE%lE",&dum_double, 

&(sll[num_points].rl),&(sll[num_points].im)); 
if(fgets(buffer,80,input_files21)!=NULL)  { 

sscanf(buffer/'%lE%lE%lE'',&freq[num_points], 

&(s21[num_points].rl),&(s21[num_points].im)); 

/^reference  plane  adjustments*/ 

sll[num_points]=cmult(sll[num_points],cexp(cmplx(0.0, 

2.0*PI*freq[num_points]*2*portl_offset/2.998e8))); 

s21[num_points]=cmult(s21[num_points],cexp(cmplx{0.0, 

2.0*PI*freq[num_points]*(portl_offset+port2_offset) 

/2.998e8))); 

s_params_to_gamma_p(d,freq[num_points],sl  1  [num_points], 
s21[num_points],&dum_complex,&(p[num_points])); 

} 

1 

} 

unwrap_phase(num_points,p,p_phase); 

fprintf(output_file/'Frequency\tmu’\tmu’'\teps’\teps"\tsigma\t 

diss_ratio\n"); 

for(loop=3;loop<=num_points-2;loop4"f)  { 

/^calculate  the  experimetal  group  delay  by  a  linear  approx  of  the  phase 
curve  over  5  points*/ 
phase_sum=0.0; 
freq_sum=0.0; 
freq2_sum=0.0; 
product_sum=0.0; 
for(loop2=-2;loop2<=2;loop2++)  { 

phase_sum+=p_phase[loop+Ioop2]; 

freq_sum4-=freq[loop+loop2]; 

freq2_sum+=freq[loop+loop2]*freq[loop+loop2]; 

product_sum+=freq[loop+loop2]*p_phase[loop+loop2]; 

} 

slope=(phase_sum*freq_sum-product_sum*5.0)/ 
(freq_sum*freq_sum-freq2_sum*5.0); 
exper_group_delay =-slope  /  (2.0*PI); 

/*loop  through  branch  cuts  to  see  which  matches  the  experimental  group 
delay*/ 
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branch=0; 

temp_sll  =sl  1  [loop]; 
temp_s21  =s21  (loop]; 

s_params_to_mateHaLparams(d,freq[loop],temp_sn,temp_s21, 

&mu_r,&epsilon_r,branch); 

branch_group_delay=group_delay(d,epsilon_r,mu_r); 

agreement= 

fabs((exper_group_delay-branch_group_delay)/ 
exper_group_delay)<  AGREEMENT; 

if(agreement)  { 

fprintf(output_file;'%E\t%7.4f\t%7.4f\t%7.4f\t%7.4f\t%E\t%E\n", 

freq[loop],mu_r.rl,mu_r>im,epsilon_r.rb-epsiIon_r.im, 

-2.0*PI*freq[loop]*8.854e'12*epsilon__r.im, 

epsilon_r.im/mu_r.im’^cmag(cdiv(mu_r,epsilot\_r))); 


)  else  { 

branch=l; 
do  { 

s_params_to_material_params(d,freq[loop],temp_sll,temp_s21, 

&mu_r,&epsilon_r,branch); 

branch_group_delay=group_deIay(d,epsilon_r,mu_r); 

agreenient= 

fabs((exper_group_delay-branch_group_delay)/ 

exper_group_deIay)<AGREEMENT; 

if(agreement)  { 

fprintf(outputJile,"%E\t%7.4f\t%7.4f\t%7.4f\t%7.4f\t%E\t%E\n", 
freq[loop],mu_r.rl,mu_r.im,epsilon_r.rV-epsilon_r.im, 
-epsilon_r.im/ epsilor  _r.rl, 
-2.0*PI*freq[loop]*8.854e-12*epsilon__r.im, 
epsilon_r.im/mu_r,im*cmag(cdiv(mu_r,epsilon_r))); 
break; 

} 

branch*=-l; 

if(branch>0)  branch ++; 
if(branch>=10)  { 

printfC'Unable  to  process  data  at  frequency=%E\n", 
freqdoop]); 
break; 

} 

1  whiledagreement); 


72 


} 

} 

stop_program: 

if(freq!=NULL)  free(freq); 
if (p_phase  !=NULL)  free(p_phase) ; 
if(sll!=NULL)  free(sll); 
if(s21!=NULL)  free(s21); 
if(p!=NULL)  free(p); 

if(input_filesl  1  !=NULL)  fcIose(input_filesl  1); 
if(input_files21!=NULL)  fclose(input_files21); 
if(output_file!=NULL)  fdose(output_file); 

} 


void  unwrap_phase(num_points,in_arr,out_arr) 
register  int  num_points; 
register  complex  '^in_arr; 
register  double  *out_arr; 

{ 

register  int  loop; 

register  double  pi,curr_phase,prev_phase,accumulate,twopi; 


pi=PI; 

twopi=2.0’^pi; 


accumulate=0.0; 
prev_phase=cpha(in_arr[l  ]); 
out_ar  r  [  1  ]  =pre  v_phase; 

for(loop=2;loop<=num_points;loop++)  { 
curr_phase=cpha(in_arr[loopl); 
if(curr_phase>=prev_phase+pi)  accumulate-=twopi; 
else  if(curr_phase<=prev_phase-pi)  accumulate+=twopi; 
out_arr[loopl=curr_phase+accumulate; 
prev_phase=curr_phase; 

} 

} 


double  group_delay(d,epsilon_r,mu_r) 
double  d; 

complex  epsilon_r,mu_r; 

{ 
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double  templ,temp2,temp3,temp4,temp5,temp6,temp7,temp8; 


tempi  =mu_r.rl*epsilon_r.rl; 
temp2=mu_r.rl*epsilon_r.im; 
temp3=mu_r.im*epsilon_r.rl; 
temp4=mu_r.im*epsilon_r.im; 

temp5=templ  -temp4; 

temp5*=temp5; 

temp6=temp2+temp3; 

temp6*=temp6; 

temp7=temp5-temp6; 

if(temp7<=0.0)  return  0.0; 

temp7=sqrt(temp7); 

temp8=templ  -temp4+temp7; 

if(temp8<=0.0)  return  0.0; 

temp8=sqrt(temp8) ; 

return  (temp8*d  /  (2.998e8*sqrt(2.0))); 


I* - V 

void  material_params_to_s_params(d,freq,sll,s21,mu_r,epsilon_r) 
double  d,freq; 

complex  *sll,*s21,mu  r,epsiIon  r; 

{ 

double  omega; 

complex  gamma,p,cl,z2,prop_const2,ctempl,ctemp2, 
ctemp3,ctemp4,ctemp5,ctemp6; 

/■^constants*/ 

cl  =cmplx(l. 0,0.0); 

/^calculated  values*/ 
omega=2.0*PI*freq; 

/*calculate  intrinsic  impedance  of  region  2*/ 
ctempl=cdiv(mu_r,epsilon_r); 
z2 =csqr  t(ctemp  1 ) ; 

/*calculate  propagation  constant  of  region  2*/ 
ctempl  =cmult(mu_r,epsiIon_r); 
ctemp2=csqrl(ctempl ) ; 

prop_const2=cmult(cmplx(0.0,omega/2.998e8),ctemp2); 
/*calculate  p*/ 
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ctempl=cmult(rmplx(-d,0.0),prop_const2); 

p=cexp(ctempl); 

/  "^calculate  gamma*/ 
ctempl =csub(z2,cl ); 
ctemp2=cadd(z2,cl); 
gamma=cdiv(ctempl,ctemp2); 

/  “^calculate  sll*/ 

ctempl =cmult(p,p); 
ctemp2=cmult(gamma,gamma); 
ctemp3=csub(cl  ^ctempl  ); 
ctemp4=cmult(ctemp3,gamma); 
ctemp5=cmult(ctemp  1  ,ctemp2) ; 
ctemp6=csub(cl  ,ctemp5) ; 

*sl  1  =cdi  v(ctemp4,ctemp6); 

/“^calculate  s21*/ 

ctempl  =cmult(gamma,gamma); 

ctemp2=cmult(p,p); 

ctemp3=csub(cl,ctempl); 

ctemp4=cmult(ctemp3,p); 

ctemp5=cmult(ctempl,ctemp2); 

ctemp6=csub(cl,ctemp5); 

*s21  =cdiv(ctemp4,ctemp6); 


} 


void  s_params_to_materiaI__params(d,freq,sl  1  ,s21  ,mu„r,epsilon_r,branch) 
double  d,freq; 

complex  si  l,s21/mu_r/epsilon_r; 
int  branch; 

complex  gamma,p,chi,x,y,cl,ctempl,ctemp2, 
ctemp3,ctemp4,ctemp5,ctemp6,ctemp7; 

/^constants*/ 

cl =cmplx(l  .0,0.0); 

/■^calculate  gamma*/ 

ctempl  =csub(cmult(sl  l,sl  1  ),cmult(s21  ,s21 )); 
ctemp2=cadd(ctempl,cl); 
ctemp3=cmult(cmplx(2.0,0.0),sll); 
chi=cdiv(ctemp2,ctemp3); 
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ctemp4=csub(cmult(chi,chi),cl); 
ctemp5=csqrt(ctemp4); 
ctemp6=cadd(chi,ctemp5); 
ctemp7=csub(chi,ctemp5) ; 
if(anag(ctemp6)'^=l)  gamma=ctemp6; 
else  gamma =ctemp7; 

/■^calculate  p*/ 

ctempl=cadd(sll,s21); 
ctemp2=csub(ctempl, gamma); 
ctemp3=cmult(ctempl, gamma); 
ctemp4=csub(cl,ctemp3); 
p=cdiv(ctemp2,ctemp4); 

/^calculate  x*/ 

ctempl=cadd(cl,gamma); 
ctemp2=csub(cl, gamma); 
x=cdiv(ctempl  ,ctemp2); 

/  “^calculate  y*/ 

ctempl  =cmplx(0.0,-2.998e8/  (2.0’^PI*freq*d)) ; 

ctemp2=cdiv(cl,p); 

ctemp3=cln(ctemp2,branch); 

y=cmult(ctempl,ctemp3); 

if(y.rl<0.0)  y=cmult(y,cmplx(-l  .0,0.0)); 

/■^calculate  mu  and  epsilon*/ 
*mu_r=cmult(x,y); 

*epsilon_r=cdiv(y,x); 

} 


void  s_params_to_gamma_p(d,freq,sll,s21,gamma,p) 
double  d,freq; 

complex  sll,s21,*gamma,*p; 

f 

double  omega; 

complex  chi,x,y,cl,ctempl,ctemp2,ctemp3,ctemp4,ctemp5,clGmp6,ctGmp7; 

/^constants*/ 

cl  =cmplx(l. 0,0.0); 

/^calculated  values*/ 
omega=2.0*PI*freq; 
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/^calculate  gamma*/ 

ctempl  =csub(cmult(sl  l,sl  1  ),cmult(s21  ,s21)); 

ctemp2=cadd(ctempl,cl ); 

ctemp3=cmult(cmplx(2.0,0.0)/Sll); 

chi=cdiv(ctemp2,ctemp3); 

ctemp4=csub(cmult(chi,chi),cl); 

ctemp5=csqrt(ctemp4); 

ctemp6=cadd(chi,ctemp5); 

ctemp7=csub(chi,ctemp5); 

if(cmag(ctemp6)<=l)  *gamma=ctemp6; 

else  *gamma=ctemp7; 

/^calculate  p*/ 

ctempl  -cadd(sl  l,s21 ) ; 

ctemp2=csub(ctempl,*gamma); 

ctemp3=cmult(ctempl,*gamma); 

ctemp4=csub(ci,ctemp3); 

*p=cdiv(ctemp2,ctemp4); 

} 


void  gamma_p_to_s_params(d,freq,gamma,p,sll,s21) 
double  d,freq; 

complex  gamma,p,*sll/s21; 

{ 

double  omega; 

complex  cl, ctempl ,ctemp2,ctemp3,ctemp4,ctemp5,ctemp6; 

/^constants*/ 

cl  =cmplx(l. 0,0.0); 

/^calculated  values*/ 
omega=2.0*PI*freq; 

/*calculate  sll*/ 

ctempl  =cmult(p,p); 

ctemp2=cmult(gamma,gamma); 

ctemp3=csub(cl, ctempl); 

ctemp4=cmult(ctemp3,gamma); 

ctemp5=cmult(ctempl,ctemp2); 

ctemp6=csub(cl,ctemp5); 

*sll=cd’  ;ctemp4,ctemp6); 


/■^calculate  s21*/ 

ctempl =cmult(gamma,gamma); 
ctemp2=cmult(p,p) ; 
ctemp3=csub(cl , ctempl ); 
ctemp4=cmult(ctemp3,p); 
ctemp5=cmult(ctempl,ctemp2); 
ctemp6=csub(cl  ,ctemp5); 

*s21  =cdiv(ctemp4,ctemp6); 


} 


void 

gamma_p_to_material_params(d,freq,gamma,p,mu_r,epsilon_r,branch) 
double  d,freq; 

complex  gamma,p,*mu_r/epsilon_r; 
int  branch; 

{ 

double  omega; 

complex  X,y,cl,ctempl,ctemp2,ctemp3,ctemp4,ctemp5,ctemp6,ctemp7; 

/^constants*/ 

cl =cmplx(l. 0,0.0); 

/■^calculated  values*/ 
omega=2.0*PI*freq; 

/■^calculate  x*/ 

ctempl  =cadd(cl, gamma); 
ctemp2=csub(cl  ,gamma) ; 
ctemp3=cdi  v(ctempl  ,ctemp2); 
x=cmult(ctemp3,ctemp3) ; 

/■^calculate  y*/ 

ctempl=cmplx(3.0e8/(omega*d),0.0); 

ctemp2=cdiv(cl,p); 

ctemp3=cln(ctemp2,branch); 

ctemp4=cmult(ctempl,ctemp3); 

ctemp5=cmult(ctemp4,ctemp4); 

y=cmult(cmplx(-1.0,0.0),ctemp5); 

/^calculate  mu*/ 
ctempl  =cmult(x,y); 

*mu_r=csqrt(ctempl); 
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/^calculate  epsilon*/ 
ctempl  =cdiv(y,x); 
*epsilon_r=csqrt{ctempl); 


} 


typedef  struct  { 
double  rl,im; 
)  coraplex; 


complex  cmplx(rl/r2) 
double  rl,r2; 

{ 

complex  temp; 

temp.rl=rl; 
temp.im=r2; 
return  temp; 

} 

double  cmag(zl) 
complex  zl; 

{ 

return  sqrt(zl  .rl*zl  .rl+zl  .im*zl  .im); 

) 

double  cmag2(zl) 
complex  zl; 

{ 

return  (zl  .rl*zl .rl+zl  .im*zl  .im); 

} 

double  cpha(zl) 
complex  zl; 

{ 

double  pha; 


if(zl.rl==0.0)  { 

if(zl.im==0.0)  return  0,0; 
if(zl.im>0.0)  return  PI/2.0; 
return  -PI/2.0; 

) 

pha=atan(zl  .im  /  zl  .r  1); 
if(zl.rl<0.0)  return  (pha+PI); 
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return  pha; 


complex  cadd(zl,z2) 
complex  zl,z2; 

{ 

complex  temp; 

temp.rl=zl  .rl+z2.rl; 
temp.im=zl  .im+z2.im; 
return  temp; 


complex  csub(zl,z2) 
complex  zl,z2; 

{ 

complex  temp; 

temp.rl=zl  .rl-z2.rl; 
temp.im=zl.im-z2.im; 
return  temp; 


complex  cmult(zl,z2) 
complex  zl,z2; 

{ 

complex  temp; 

temp.rl=zl.rl*z2.rl-zl.im*z2.im; 
temp.im=zl  .rl*z2.im+zl  .im*z2.rl; 
return  temp; 


complex  cdiv(zl,z2) 
complex  zl^2; 

{ 

complex  temp; 
double  rtemp; 

rtemp=cmag2(z2); 

temp.rl={zl  .rl*z2.rl+zl  ,im*z2.im)  /  rtemp; 
temp. im=(zl  .im*z2  .rl-zl  .rPz2.im)  /  rtemp; 
return  temp; 


complex  cexp(zl) 
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complex  zl; 


{ 


} 


complex  temp; 
double  decay; 

decay =exp(zl.rl); 
temp.rl=decay  *cos(zl  .im); 
temp.im=decay*sin(zl  -im); 
return  temp; 


complex  conj(zl) 
complex  zl; 

{ 

complex  temp; 


] 


temp.rl=zl.rl; 
temp.im=-zl.im; 
return  temp; 


complex  cp2r(zl) 
complex  zl; 

{ 

/♦when  sent  a  polar  complex  with  magnitude  in  real  part  and 
phase  (radians)  in  imaginary  part,  this  returns  a  rectangular  complex*/ 

complex  temp; 

temp.d=zl  .rl*cos(zl  .im); 
temp.im=zl  .rl*sm(zl  .im); 
return  temp; 

} 


complex  cr2p(zl) 
complex  zl; 

{ 

/♦when  sent  a  rectangular  complex,  this  returns  a  polar  complex  with 
magnitude  in  real  part  and  phase  (radians)  in  imaginary  part*/ 

complex  temp; 

temp.rl=cmag(zl  ); 
temp.im=cpha(zl); 
return  temp; 
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complex  csqrt(zl) 
complex  zl; 

{ 

complex  temp; 


temp=cr2p(zl); 
temp.rl=sqrt{temp.rl); 
temp.im=temp.im/2.0; 
temp=cp2r(temp) ; 
return  temp; 


complex  croot(zl,n,i) 
complex  zl; 
int  n,i; 

{ 

/^returns  the  ith  value  of  zl*’^(l/n) 
i  is  between  0  and  n-1 

V 

complex  ctemp; 

register  double  mag,phase,dn; 


dn=(double)n; 

ctemp=cr2p(zl); 

mag=pow(ctemp.rl,l  .0  /  dn); 

phase=(ctemp.im+2.0*PI*i)  /  dn; 

ctemp=cmplx(mag*cos(phase),mag*sin(phase)); 

return  ctemp; 


complex  cln(zl,branch) 
complex  zl; 
int  branch; 

{ 

complex  temp; 

temp.rl=log(cmag(zl)); 
temp.im=cpha(zl)+2.0*PI*branch; 
return  temp; 


GOVERNMENT  PfllNTfNG  OFFICE 


82 


MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 

O 

Force  Command,  Control,  Communications  and  Intelligence  (C  l)  activities 
for  all  Air  Force  platforms,  /t  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  BSD  Program  Offices  (POs)  and  other 
BSD  elements  to  perform  effective  acquisition  of  C^I  systems.  In  addition, 
Rome  Laboratorys  technology  sports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


